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ABSTRACT 



We present a thermo-kinetic description of anomalous diffusion of single-particles and 
clusters in a viscoelastic medium in terms of a non-Markovian diffusion equation involving 
memory functions. The scaling behaviour of these functions is analyzed by considering 
hydrodynamics and cluster-size space random walk arguments. We explain experimental 
results on diffusion of Brownian particles in the cytoskeleton, in cluster-cluster aggregation 
and in a suspension of micelles. 
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I. INTRODUCTION 



Diffusion is an irreversible phenomenon which arises in systems of different natures. The 
understanding of the mechanisms intrinsic to this process remains an open question in cases 



where the host liquid displays a variety of time and length sea 



in slow relaxation systems 



soft matter aggregations 



9 



3. 



10 



4 [jin viscoelastic media 



11 



es. Such cases include diffusion 



6|, Lfl, ISj or in the presence of 



121 ] . To find out the peculiar behaviour of the mean 
square displacement (MSD) of the particles observed in these situations is a theme of vital 
importance in areas such as soft condensed matter, biophysics and material science. 

The properties of viscoelastic fluids can be characterized through microrheological exper- 



iments in which the MSD of passively diffusing particles is measured 



it has been shown both experimentally 



4, la, [3]. For example, 



that the MSD of a probe 



4| and theoretically 

particle moving through a medium made up of a solution of F-actin networks obeys a power 
law, with an exponent depending on the aspect ratio between the particle radius and the 
characteristic length of the polymer network. Similar results for the MSD have also been 

n 

obtained in the case of living cells [6] and then used to explain the viscoelastic properties 
of the intracellular medium. This explanation is usually given in terms of the single-point 



generalized Stokes-Einstein relation in the frequency domain 



13 



14 



1 51 ] . The viscoelastic 



and non-Newtonian properties are modified when the system is subjected to shear stresses 



or temperature gradients 



16( , and to geometrical confinement 



171 ] in which the finite size 



of the particles may also play an important role in the dynamics of the system 



Matter aggregation phenomena in viscoelastic environments show up many other interest- 



ing 



id, 



eatures both at molecular as well as supramolecular or mesoscopic levels of organization 



19 



201 ]. Although the molecular level is by now thoroughly accessible through modern 



4 



experimental techniques and computer simulations, its mesoscopic counterpart lends itself 
readily amenable to theoretical modeling, regardless of whether the process under study re- 
flects subdiffusive or superdiffusive dynamic properties 



mention formations of networks and gels 



id, 



8 



15| . As for the former, one could 



2l|, viscoelastic phase separation 



12 



well as formations of biomembranes and/or soft bio materials 



22| as 



9] . As for the latter, one could 



point out colloidal and protein systems in which nucleation-growth phase transformations 



often lead to a global superdiffusive behavior, especia! 



tein crystals or colloidal crystals comes into play 



18 



ly when formation of non-Kossel pro 



23] . In such systems, finiteness 



id 



and mass flucutations [24j of the basic constituent (a large Brownian-type particle, or a clus- 
ter of smaller particles) as well as some detected viscoelastic response of the system upon 
aggregation are of prior significance and point undoubtedly to more realistic descriptions. 
In this paper, we propose a thermo-kinetic framework based on mesoscopic nonequilib- 



rium thermodynamics 19J giving a coarse-grained description of the anomalous transport 
of particles and their clusters in a viscoelastic medium through a non-Markovian diffusion 
equation, engaging thoroughly the slowest dynamic variable termed the reaction coordinate 
20t ] . This equation is derived under the assumption of a non- instantaneous response of the 



25 



26 



271 ] . These coeffi- 



system in which transport coefficients become memory functions 
cients are introduced into the description throu gh t he expression of the probability current 



19 



20]. 



obtained from the entropy production equation 

The paper is organized as follows. In section II we will derive the non-Markovian diffusion 
equation accounting for the dynamics of the particles in a viscoelastic medium. In Sec. Ill 
we will analyze the scaling of the memory function in a polymer solution and in Sec. IV 
we will analyze the scaling of the MSD in cluster-cluster aggregation. Section V will be 



devoted to comparing theory with experiments and Sec. VI to discussion and summary of 



5 



our main results. 



II. THE NON-MARKOVIAN DIFFUSION EQUATION 



Let us consider the Brownian motion of a particle through a viscoelastic medium such as a 
polymer melt in which the macromolecules are distributed as to form a network which exerts 
elastic forces on the particle 5|, la, LZI]. In general, the medium may also involve dynamical 
processes associated with changes in the conformation and the growth of its costitutents 
as is the case of the polymerization of microtubules or the aggregation and crystallization 



of lipids and proteins 



10 



231 ]. For simplicity we will obtain the properties of this 



viscoelastic heat bath under an effective medium a ppr oximation 



19 



20 



by means of the probability 



The dynamics of the particle can be described 
density p(r, t) depending on its position r and time t. Since p is normalized to the unity it 
obeys the continuity equation 



0_ 

Of 



P 



-VJ, 



(1) 



where V = . The explicit form of the diffusion current J can be inferred from the entropy 
production of the diffusion process. This quantity can be obtained from the nonequilibrium 



entropy |19 



20. 



28| 



sit) = -k B / pir, t) In — dr + s t 

J Pie 



(2) 



where ks is Boltzmann's constant and si e its value at the reference local equilibrium state 
characterized by the probability density 



(3) 



Here /z/ e is the corresponding chemical potential and 0r(r) the interaction potential which 
contains two contributions: <p e related to the potential energy of the external forces and 



b (r) resulting from the interactions between the particle and the viscoelastic medium. 

By taking the time derivative of Eq. (121), using (CEJ) and integrating by parts we obtain 
the entropy production 

^/j-V/idr, (4) 



°- T 



where we have defined the nonequilibrium chemical potential: /i(r, t) = k B T In p/ pi e + pi e . A 
similar type of chemical potential was obtained during a structural transition of mono-olein 



29] 



bicontinuous cubic phase induced by inclusion of protein lysozyme solutions, cf. Ref. 
In this case, the ratio of nonequilibrium and equilibrium densities pj p\ e is interpreted as a 
confinement factor which corresponds to the ratio of the volume of water channels, v w , and 
the volume available for the center of lysozyme molecules v a , namely v w /v a , cf. Table II 
and Eq. (9) therein. The structural frustration is caused by a moderate attraction of the 



lysozyme molecules mediated by the lipid system 



phase separation 



121. 



9j, and can be attributed to a (viscoelastic) 



we may establish a 



Following now the scheme of nonequilibrium thermodynamics 
relationship between the probability current J and its conjugated force Vp. Taking into 
account that in a viscoelastic medium the response of the system is, in general, not instan- 



taneous 



251 ] we obtain 

J(r, t) = - f C(t - t')p(r, f ) (V0 r + k B TV lnp(r, f)) dt', (5) 
Jo 

where we have used the explicit expression of p\ eq and introduced the memory function 
C(t — t'). Notice that in the case of non-isotropic systems £(£ — t') becomes a tensor. 
After substituting Eq. (jSJ) into ([T]), one obtains the non-Markovian diffusion equation 

= f C(t - V • [p(r, V0 T ] dt' + f D(t - t')V 2 p(r, t')dt', (6) 
at Jo Jo 



where we have defined the memory function D(t) = /cgT£(£). The first term on the right- 
hand side of Eq. ([6]) contains the force exerted by the viscoelastic medium on the particle. 
In the case of diffusion in an intracellular medium, there may be elastic forces arising from 



the polymer network and the forces due to the 



activity of molecular motors 



18 



30|. In 



the case of diffusion-limited growth processes 
nature or be related to the presence of electrostatic fields 18 
obeys the Stokes-Einstein law 



9j, these forces could also be of an elastic 



231 ] . The diffusion coefficient 



(7) 



where (3{u) = 6nari(uj) is the friction coefficient containing the radius a of the particle and 



the frequency-dependent viscosity rj(w), 



15 



3l|. The memory function and (3(uj) are 
related through the inverse Laplace transform: ((t) = £~ l [/3 -1 (u;)]. Equation is similar 
to the ones previously derived in the context of the continuous-time random walk formalism 



and of projector operator techniques 



25 



27] - 



Explicit expressions for the memory functions ((t) and D(t) can be obtained from hydro- 



dynamics. It has been shown 



331 ] that the frequency- dependent corrections to the friction 
coefficient are given by a power expansion of the inverse penetration length: ao = (lu/uq) 1 ^ 2 , 



with vq the kinematic viscosity of the heat bath 



j|. In the case when the heat bath is 



a viscoelastic fluid, the inverse penetration length incorporates the frequ ency 



kinematic viscosity v(uS): a(u) = [to / 'v(uj)] 1 I 2 , leading to the expression 



/3(u) ~ A) [1 + aa(u)] 



14 



32 



dependent 



33] 



(8) 



This length behaves as: a(u>) ~ (aA w ) 1 (t d oj)( 2 ), where td is a characteristic re- 
laxation time. We have assumed that, within a given range of frequencies, the kinematic 



viscosity of the viscoelastic medium obeys the scaling: v(oj) ~ u (tdu) 



14J. The quanti- 



ties S and A w may, in general, depend on the free space around the particle [34| and on the 



temperature 



14 



35]. The scaling factor has its origin in the fact that the relaxation 



time of the system may depend on the concentration of, for instance, polymers or micelles 



making up the viscoeastic medium 



14 



35 



36]. 



The exponent 5 also admits different interpretations. An important one of these is the 
one related to the Brownian motion of clusters or small crystals which are still in their late 



growing stage, e.g. when forming semi-orderly amorphous mesostructures 



34J. In this case, 



the exponent may reflect the specific dynamics obeyed by the radius of these objects as has 



been analyzed in Ref. 



101. 



III. SCALING OF THE MEMORY FUNCTION 



The scaling of the memory function and, consequently, the microrheological properties 
of the viscoelastic medium can be inferred from the MSD of a test particle. The MSD can 
be obtained experimentally by using diffusing-wave spectrometry techniques or video-based 



methods 



6. 



13]. It has been found that, in a certain range of frequencies, the complex 



shear modulus obeys the scaling behavior: G'^J^u) ~ uj 1 with 7 < 1, 

These experimental results can be described by using Eq. (jH]) in the case of 4>t — 0. For 
frequencies lying in the range r^ 1 < u <C /3q, the leading term of the Eq. ([8]) yields 



The inverse Laplace transform of (Q can be performed by using the Tauberian theorem 
One then obtains the memory function 



(9) 



22]. 



9 



where T(b) is Euler's gamma function. After substituting (1101) into ([H]) with <$t = and 
using the inverse Laplace transform of Eq. (J7J), one arrives at the non-Markovian diffusion 



equation 



dp 

~dt ~ tdT 



Dn 



1-8- 



t t-t' 



■ 1+6 - 

, 2 . 



\ T D 



V 2 p(r,t')dt', 



(11) 



27]. 



which is similar to the one previously derived in Ref. 

It is important to realize that equation (fTTj) can be recast in the form of a fractional 
diffusion equation 37] by Laplace transforming it and multiplying the result by (std) -1 - 
Then, taking the time derivative of the inverse Laplace transform of the resulting expression, 
one obtains 

^(r,*)= oA 1- *A>V a p(M), (12) 
in which the exponent is a = (3 — <5)/2, with 1 < 5 < 3 and the integro-differential operator 

,D}-*A(t') = T-\a)- / (-*-) A(t')dt' (13) 



dtJo \t-t 



371 ] . This feature demon- 



entering Eq. (fT2l is the zero order Riemann-Liouville operator 
strates that our method based on a scaling of the memory function and a description using 
fractional derivatives leads to an identical kinetic equation. 



The value of the exponent 5 characterizing the properties of t 
in general be obtained by analyzing the hydrodynamic effects 



re viscoelastic medium can 
For simplicity, we will 



show how this exponent may be inferred from the evolution in time of the MSD 

(r 2 )(t) = J r 2 p(r,t)dr. (14) 

The evolution equation for (r 2 )(t) follows by taking the time derivative of Eq. f[T4T) and then 
using Eq. ffTTj) . After integrating by parts, one obtains 



dt 



D 



t \ — 

Td - 



(15) 



10 



whose solution is 



(r 2 )(t) ~ t~ 



(16) 



Equation ffTBT) can now be compared with the experimental result (r 2 ) exp ~ ti, obtained 
for a Brownian particle in the intracellular medium |6j, yielding 5 = 3/2. In view of Eq. 
(jHJ), the mobility is given by (3(oj) — uj~±. Since (3(oj) oc rj e ff(u) with rj e ff(u) the effective 
viscosity of the medium and G" msd {oj) ~ ^Veff, it follows that our result is in good agreement 
with the experimental result G" MSD ~ u*, reported in Ref. [13]. 



IV. SCALING OF THE MSD IN CLUSTER-CLUSTER AGGREGATION 



Cluster-cluster aggregation in a viscoelastic medium at high temperatures can be viewed 



as a (^-dimensional random walk (RW) in the space of cluster sizes x, 10j. Below, we present 
some arguments which lead to the recovery of the experimentally justified result given by 
Eq. (flfil) also at the level of clusters' RW. 

A coupling of late-time growing and mechanical relaxation conditions has been proposed 
in Ref. [10j in order to properly describe the clusters' RW in the fluctuating viscoelastic 



medium. The role of MSD for such RW is played by t 



squared, R av 2 (t), which scales with the time t as 



re average radius of the clusters but 



10 



2l| 



iC 2 (*)~* 2/(,m) . 



(17) 



According to Ref. lOj], and based on the above scaling relation, for a planar geometry 
the value of the exponent of the MSD is vmsd = 2/3, which is the harmonic mean of two 
other exponents obtained for d — 1 and d = 3, namely vmsd = 1 and Umsd = 1/2, respec- 
tively. This is due to some physical interpretation that the harmonic mean well describes 



an optimal (excluded-volume involving 



34|) kinetic pathway for a clustering thereby intro- 



11 



duced, which ultimately favors d = 2. This interpretation is in accordance with soft-glass 
matter-organization conditions which favor characteristic exponents less than 3/4, see Ref. 



3l| and references therein. If we, however, take an arithmetic instead of the harmonic mean, 
and apply it to both extreme exponents, umsd = 1 and v MSD = 1/2, which looks like a 



customary fitting procedure applied by experimenters, cf. Ref. 



1[, we finally arrive at the 



arithmetic-mean exponent: vmsd = 3/4, which is characteristic of the viscoelastic matter 



relaxation (under late-stage growing conditions) 
view offered by the present study. 



101 ] . This relation supports the general 



The scaling formula (ITT)) can be used to obtain the particle MSD of specific concentrated 



colloida 
in Ref. 



systems. For example, in the case of the sol-gel continuous phase transition reported 
the exponent of the initial and final phases can be acce ptably well- reproduced by 



id, 



ii 



2JJ . In our model, the 



taking into account the pivotal role of clusters' surface a^ -1 )/^ 
characteristic exponent v c h ar of these phases correspond to v c h ar — 1 for d — 1 (no surface) 
and Vchar = 0.66 for d = 2 (a planar geometry) which is very close to the measured 0.7, see 



Figure 4 of Ref. 



Thereby, experiment and theory show almost the same exponents 



values under which the same dynamic behavior, i.e. a passage from a diffusive (sol) to some 



subdiffusive (gel) state can be seen 



cf. Figure 3 of the review 



lOj ] . This is also confirmed by relaxation measurements, 



and our considerations below. 



The MSD R 2 av introduced in Eq. (ITTj) can be further related with the creep compliance ht(t), 
defined by the average strain to be assigned to a probe particle immersed in the viscoelastic 



medium (e.g., the cytoskeleton) , divided by the corresponding shear stress [31 



36] 



R 2 av (t) OC K(t), 



In the microrheological high-temperature environment considered in Ref. 



(18) 



10], we have 



made an attempt to define in a phenomenological way the inverse of k, namely 



12 

the internal matrix stress cxm(£) accumulated within the inter-cluster spaces of the matrix, 
assuming that in the late-stage growing conditions the average strain can be set constant. 
Taking this relation to be valid, the creep compliance for our high T clustering viscoelastic 
system satisfies n(t) ~ a M ~ 2 (t), which after taking into account Eqs. (fT7|) and ffl8l) leads to 

a M (t)oct-^ d+1 \ (19) 

i.e. a slightly modified result than while entirely based on the before applied phenomenology 



lOt ] . From this relation, it follows that the scaling of the internal matrix stress will, in 
general, be a function of the relevant degree of freedom, d -the space dimension, in the same 
way as the MSD. 

V. COMPARISON WITH EXPERIMENTS 

A. The cytoskeleton response: fast-slow crossover of the complex shear modulus 

Recent experiments measuring the viscoelastic response of the cytoskeleton through its 
complex shear modulus G"(oj) have shown a crossover from the scaling behaviour G"{uj) 



!]• 



^3/4 f or } ar g e values of the frequency to an almost constant value for low frequencies 
These experiments were performed by attaching a magnetic bead to the cytoskeleton network 



!]• 



and then applying a periodic magnetic field that exerts a torque on the particle 

This crossover can be explained by means of the theory presented in section II if one 
assumes that the total force exerted on the attached probe particle can be modelled by 
means of the harmonic force F(r) = — uj^t with ujq a characteristic frequency, and that in 
Eq. (JTUD 5 = 3/2. 

From our description one may show that the complex shear modulus is related with 



13 



(f 2 )(u;), through the well-known expression 31 1 



G"(u) 



k B T 1 

na uj (r 2 ) 



(20) 



where we have used the relation r)(ou) = r/oC obtained after taking the Laplace transform 
of Eq. ([6]) and using ([7j). From the resulting equation we can obtain the MSD 



<r 2 >M 



k B T 1 



1/4, ,3/4 



Td w 



(21) 



Substituting now this expression into Eq. (1201) . one finds the complex shear modulus of a 
probe particle attached to the network 



G"(u) = n (2^ + r^V/ 4 ) 



(22) 



As shown in Figure 1, our result (1221) is in good agreement with the experiments reported 



in Ref. 



1[ . The linear model introduced in this subsection proves to be adequate in order to 



reproduce experiments performed at constant temperature. The effect of temperature can 



be incorporated into the description 



librium distri 
time td 



14 



Dution function (j3J) 
331. 



16 



w 



taking into account entropic effects in the local equi- 



171 ] and through the dependence on T of the relaxation 



B. Brownian motion in a solution of giant micelles 



We will now proceed to analyze the behavior of the MSD of a particle undergoing Brown 



ian motion in a semidilute solution of worm-like micelles 



39] . The host fluid will be modelled 



as a structured viscoelastic medium acting on the particle with the non-linear force 



F e \ x ) = -V0 = F cos(A _i a; 



(23) 



14 



where F is the amplitude of the force, considering the one-dimensional case for simplicity. 
A can be interpreted as the characteritic distance between the particle and the cage formed 
by the micelles. 



391 ] . we will ap- 



In order to describe the experimental results for the MSD reported in 
proximate p(r, t + z) through: p(r,t+z) ~ p(r, t) + (dp/dt)(t)z + 0(z 2 ). This approximation 
takes into account that for the subdiffusive processes we are considering here, the relaxation 
of the system is not exponential and thus no clear separation of times scales can be assumed. 
As a consequence, at larger times the dynamics of the system remains non-Markovian, but 



with time-dependent coefficients 



271 ] . Thus, keeping the lowest order term of the expansion 



and substituting the resulting expression into ([6]), one obtains 



d_ 

Of 



q2 

p{x,t) = D{t)—p{x,t) + ({t)F 



d_ 

dx 



p(x, t) cos (A 1 x^j 



(24) 



where we have introduced the time-dependent coefficients 



((f) = / C(z)dz and D{t) = k B T£(t), 
Jo 



(25) 



where the memory function ((t) is given by Eq. (flUj) . From Eq. (1241) we obtain the evolution 
equation for the MSD 



,i (x 2 )=2D -2^(xcos(A- 1 x)), 



dr 



(26) 



where we have introduced the new variable r(t) = / ((t)dt. This equation implies that the 
effects due to the non-linear elastic force on the behaviour of the MSD appear through the 



171 ] . Since the solution of 



average work (xF el ) which this force exerts over the particle 
Eq. (124"]) with F el given by ff23]) cannot be obtained analytically, the average (xcos(x/X)) 
cannot be explicitly computed. However, we will assume small deviations with respect to 



the average value of the position of the particle, and then use the following expansion around 



15 



the mean value (x)(t): (f (x)) ~ / ((x)) + O ((x 2 — (x) 2 )). The resulting equation is 



-^(x 2 )(r) = 2D - 2j^(x)(r) cos [A-^Xt) 



which must be solved together with the approximated evolution equation for (x)(r) 



(27) 



FqDq 

' kf/T 



cos 



(28) 



Integrating Eqs. ( 12T1) and ( 1281) one can obtain the following expression for the MSD 



(x 2 )(t) ~ (x 2 )(t ) + 2D r(t) +4A 2 tanh 



FqDq 

2Xk B T 



r(t) 



(29) 



where t is a cut-off time, r(t) = ^ 3 _ s ^ s _^ T -1 [(5 — l)/2] 1 )/ 2 ^(3-5)/2 an j^ f or conven j ence; 
we have taken the first term of the expansion of arctan(x) 2 ~ x 2 with x = tanh[^^r (t)) . 
The behaviour of the MSD fT29|) as a function of time (solid line ) is shown in Figure 2, 

39j. We have used the 



hat 



39| 



and compared with experimental data (circles) taken from Ref. 
following values for the parameters: Dq ~ 3 • 10 _5 /im 2 s _1 , fcgT ~ 4.2 • 10~ 6 gr fim 2 s~ 2 , 
Fo ~ 0.5grfim s~ 2 , td — 0.1s and A ~ 8 • 10 _3 /im. From Figure 2 we may conclude t 
the expression for the MSD we have obtained agrees with microrheological experiments 
within the time interval represented. Notice that the plateau is a signature of the existence 
of cage effects, which in our model become manifest at the maximum value of the elastic 
force F el . The exponent characterizing the memory function of the solution of giant micelles 
is 5 = 9/5, which differs from that of our previous example. 



VI. DISCUSSION 



In this article, we have proposed a thermo-kinetic approach to analyze anomalous diffusion 
of particles and clusters in systems having a viscoelastic response. A non-Markovian diffusion 



16 



equation has been derived from the entropy production at the mesoscale related to the 
diffusion process. 

Using hydrodynamic arguments, we have obtained an expression for the memory function 
which obeys a power law in frequencies and leads to the scaling behaviour of the mean square 
displacement of the Brownian particle. 

Cluster-cluster aggregation in the viscoelastic medium considered as a d-dimensional 



9 



23j . The role of MSD 



random walk in the space of cluster sizes x has also been analyzed 
for such a RW, played by the average radius of the clusters but squared, has been used to 
argue that the time dependence of the matrix stress and the creep compliance of the colloid- 
type system could scale with a law depending on the relevant degree of freedom, d -the 
space dimension of the Brownian particle. Moreover, we have shown that the characteristic 
exponent 3/4 of the MSD can be obtained by applying an arithmetic mean to both extreme 
exponents (d — 1 and d — 3). This procedure coincides well with the one followed when 



3l|. 



obtaining the exponents from experimental data 

As particular cases of our general formalism, we have studied the dynamics of a particle 
attached to the cytoskeleton and of a Brownian particle moving through a semidilute solution 
of giant worm-like micelles. In the first case, we have explained the crossover of the complex 
shear modulus of the cytoskeleton observed in experiments in terms of an harmonic-force 
model. In the second case, we have analyzed the dynamics of a Brownian particle moving 
through a semidilute solution of giant worm-like micelles by proposing a nonlinear force 
model that takes cage effects into account. The good agreement between experiment and 
theory, shown through Figures 1 and 2, allows one to conclude that the theoretical approach 
we have proposed in this article can explain the behaviour observed in microrheological 



experiments 



31 



39]. 



17 



Our study suggests for the possibility of providing justifications of experimental results 



8 



23|. The 



on diffusion in a viscoelastic medium in terms of a Smoluchowski description 
theory proposed may constitute a useful tool ready for examining the dynamics of the 
soft-matter systems at a mesoscale. 
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Figure 1. 

Complex shear modulus of the cytoskeleton as a function of frequency for the case of 
a magnetic bead attached to the cytoskeleton network. An oscillatory magnetic field was 
applied in order to measure the response of the material. Circles represent experimental 



data taken from Ref. 



l| whereas the solid line a fit of the data obtained from Eq. (1221) . for 



a = 2.25-l(r 6 m, (3 ~ MOV 1 , r/ ~ l-l^kgrn^s' 1 , uj ~ 1.77-KrV- 1 and r D ~ 4-l(T 5 s. 



Figure 2. 



The MSD of a polystyrene particle in a semidilute solution of worm-like micelles as a 
function of time (1231) . The black circles are experimental data obtained by means of diffusing 
wave spectroscopy techniques, and taken from Ref. 39]. The solid line was obtained with 



Eq. fl29l . The values of the parameters are given in the text. 
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